Mon. Not. R. Astron. Soc. 000, [|-?? (2002) Printed 1 February 2008 (MN IAT^X style file vl.4) 



On the interaction of FR-II radio sources with the 
intracluster medium 



P. Alexander 

Astrophysics Group, Cavendish Laboratory, Madingley Road, Cambridge CBS 0HE 



Accepted ?????? . Received ??????; in original form ?????? 



(N 
O 
O 
(N 

Sh 

Oh 
< 



> 

m 

O 
(N 
O 

6 

CO 

> ! 1 INTRODUCTION 



ABSTRACT 

The effect of the expansion of powerful FR-II radio sources into a cluster environment is 
discussed. The analysis considers both the thermal and temporal evolution of the ICM which 
has passed through the bow shock of the radio source and the effect of this swept-up gas 
on the dynamics of the radio source itself. The final state of the swept-up ICM is critically 
dependent on the thermal conductivity of the gas. If the gas behind the bow shock expands 
adiabatically, and the source is expanding into a steeply falling atmosphere, then a narrow 
dense layer will form as the radio source lifts gas out of the cluster potential. This layer has 
a cooling time very much less than that of the gas just ahead of the radio source. This effect 
does not occur if the thermal conductivity of the gas is high, or if the cluster atmosphere is 
shallow. The swept-up gas also affects the dynamics of the radio source especially as it slows 
towards sub-sonic expansion. The preferential accumulation of the swept-up gas to the sides 
of the cocoon leads to the aspect ratio of the source increasing. Eventually the contact surface 
must become Ray leigh- Taylor unstable leading both to inflow of the swept-up ICM into the 
cavity created by the cocoon, but also substantial mixing of the cooler denser swept-up gas 
with the ambient ICM thereby creating a multi-phase ICM. The radio source is likely to have 
a marked effect on the cluster on timescales long compared to the age of the source. 

Key words: hydrodynamics - galaxies: active - galaxies: jets - X-rays: galaxies : 
clusters - galaxies: cooling flows 



Observations with both the ROSAT and CHANDRA satel- 
lites have demonstrated clearly that the presence of a 
powerful FR-II radio source very significantly affects the 
ICM/IGM in the vicinity of the source. The simplest model 
for such a radio source has an over-pressured cocoon of jet 
material expanding supersonically into the external gas; ex- 
panding ahead of the cocoon a bow shock sweeps up the 
surrounding gas which then forms a shell between the bow- 
shock and the contact surface separating shocked gas and 
cocoon material. In this model the cocoon should there- 
fore be devoid of ICM/IGM and the strong shock should 
heat the ICM/IGM. Some observations, but not all, sup- 
port this general picture. Both ROSAT and CHANDRA 
data for Cygnus A (Carilli, Perley & Harris 1994; Smith et 
al. 2002) show cavities in the X-ray emitting gas coincident 
with the radio cocoon and the effects of the bow shock can 
be clearly seen. The observations of both the A426/3C84 
interaction (Perseus; Fabian et al. 2000; Fabian et al. 2002) 
and A2052/3C317 show clear evidence for a radio cocoon 
devoid of X-ray emitting gas although in both cases the ra- 
dio structure is complex (Rizza et al. 2000; Blanton et al. 
2001). In two cluster/radio-source interactions in A133 and 



A2626 studied by Rizza et al. (2000) using ROSAT the X- 
ray emitting gas may well occupy the same volume as the 
radio-emitting plasma. The interaction of FR-I or complex 
radio sources with their environment has often been inter- 
preted in terms of the buoyancy effects of the low-density 
radio-emitting plasma as in M87 (Bohringer et al. 1995; 
Briigen and Kaiser 2001), Hydra A (McNamara et al. 2000a; 
David et al. 2001; Nulsen et al. 2002) and Centaurus A (Sax- 
ton, Sutherland & Bicknell 2002). A recent review of the X- 
ray observational data on radio-source/cluster interaction 
is given in McNamara et al. (2002b). This interaction be- 
tween radio source and the ICM has led a number of authors 
to consider whether radio sources could solve the so-called 
"cooling-flow problem" (e.g. Binney and Tabor 1995; McNa- 
mara et al. 2000b; Reynolds, Heinz & Begelman 2001; Chu- 
razov et al. 2001; Quilis, Bower & Balogh 2001; Bohringer 
et al. 2002), in which the cool gas expected to accumulate 
at the cluster centre as a result of cooling is not detected. 

A number of authors have considered various aspects of 
the physics and evolution of the gas swept-up by the passage 
of the radio source. Scheuer (1974) pointed out the poten- 
tial for the swept-up gas to become Rayleigh Taylor unstable 
given that the density of the cocoon must be very much less 
than that of the external gas. Gull & Northover (1973) had 
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already considered the evolution of a light bubble of gas in a 
cluster potential as a possible model for radio source confine- 
ment, while Smith et al. (1983) determined the conditions 
under which a bubble, rather than a jet, would be produced 
in the Blandford and Rees (1974) twin-exhaust model. This 
picture has recently been updated by Churazov et al. (2001), 
Briiggen & Kaiser (2001), Quillis, Bower & Balogh (2001) 
and Saxton, Sutherland & Bicknell (2002) who have investi- 
gated the evolution of buoyant bubbles. These studies show 
that the buoyant material can give rise to a mushroom-type 
structure which although susceptible to Rayleigh- Taylor in- 
stability can significantly affect a cooling flow within the 
cluster. Heinz, Reynolds & Begelman (1998; hereafter HRB), 
calculated the expected X-ray emission from the swept-up 
ICM using a simplified geometry for the radio source and as- 
suming uniform conditions between the bow-shock and the 
contact surface. A more detailed calculation of the heat- 
ing effects of a radio source on its environment was made 
by Kaiser & Alexander (1999; hereafter KA99) assuming 
an elliptical source shape, self-similar evolution and adia- 
batic conditions in the swept-up gas. They predict signifi- 
cant structure in the X-ray emission behind the bow shock 
and flow in this region away from the hotspot. Reynolds, 
Heinz & Begelman (2001) present detailed simulations of a 
radio source expanding into an atmosphere; for sources ex- 
panding supersonically their results are in agreement with 
the analytical models of KA99, but they are able to fol- 
low the evolution of the radio source into a phase in which 
it reaches pressure balance with the surrounding gas and 
again find significant structure in the expected X-ray emis- 
sion. Despite this literature a number of important problems 
concerning the evolution of the swept-up gas have not been 
fully addressed. In particular what is the long term effect 
of the radio source on its environment on time-scales long 
compared to the source life time, and what effect does the 
swept-up gas have on the evolution of the radio source. Heat- 
ing of the ICM/IGM by radio sources may have important 
and measurable cosmological effects (e.g. Yamada & Fujita 
2001). 

In this paper I extend this earlier work to determine 
the effects of the expansion of a powerful FR-II radio source 
into a cluster- like environment. The structure of this paper 
is as follows. In Section 2 I outline the main physical con- 
siderations appropriate to the analysis of the evolution of 
gas swept-up by the radio source. In Section 3 I consider the 
evolution of this gas in detail and in Section 4 I consider 
the effects of the gas on the dynamics of the radio source. 
Finally, in Section 5, I discuss the implications of these re- 
sults. 



2 GENERAL CONSIDERATIONS 

Before discussing more detailed models, it will be useful to 
consider the main factors which are likely to be important 
for the swept-up gas and the radio source itself. I begin 
by reviewing the likely physical conditions in the swept-up 
gas. Figure 1 shows a sketch of the assumed geometry of 
the source. For simplicity I shall refer to the atmosphere 
into which the radio source propagates as the Intra Cluster 
Medium (ICM). 

As the radio source expands it must do work on the sur- 



rounding medium. For supersonic expansion the ICM will be 
heated by the bow shock and is also lifted out of the poten- 
tial well in which the host galaxy sits; this swept-up gas re- 
sides between the bow shock and contact surface. The sound 
speed in this shocked gas must be of order the velocity of the 
bow shock itself, and therefore the sound crossing time must 
be somewhat less than the age of the radio source suggest- 
ing that the pressure between the bow shock and contact 
surface should be approximately constant. This argument is 
confirmed by numerical simulations of FR-II radio sources 
(e.g. Reynolds et al. 2001) and the calculations of KA99. 
There must be pressure variations within the shocked gas 
around the periphery of the source since the pressure at the 
hotspot exceeds that within the cocoon; these pressure vari- 
ations drive a flow within the swept-up gas sending material 
away from the hotspot in the direction of the host galaxy 
(KA99, Reynolds et al. 2001). 

As the ICM passes through the bow shock it is of 
course heated. As this gas flows downstream of the bow 
shock it will evolve adiabatically unless there is significant 
radiative heating / cooling or the thermal conductivity is 
high. The issue of the cooling of the gas will be consid- 
ered in detail below. The X-ray gas is optically thin, there- 
fore it is only necessary to consider the thermal conductiv- 
ity as a heating source. The thermal conductivity can be 
written in the form 1.3r]n t: X e (kT e /m e ) 1 ^ 2 , where the sub- 
script e refers to electron values, \ e is the electron mean 
free path, and the dimensionless factor r\ allows for sup- 
pression of the thermal conductivity below the canonical 
form given by Cowie & McKee (1977). The heating time, 
th, for a temperature variation of AT to be removed in a 
layer of gas of thickness Ar is given by th ~ 1.5nkATAr /q, 
where the heat flux, q = 1.3r/n e \ e {kT e /m e ) 1/2 A(kT)/ Ar. 
This simplifies since (fcT e /m e ) 1/ ' 2 ~ v e ~ c s (m p /m e ) 1//2 , 
hence th ~ 0.03i s (Ar/A e )/r/, where t s is the sound crossing 
time across the region of size Ar. The electron mean free 
path is given approximately by: 

Ae " 23 (i^) ( io-^-O lkpc - 

For Cygnus A at a radius of approximately 100 kpc n e ~ 
7 x 10~ 3 cm~ 3 , T c ~ T gas ~ 7 x 10 7 K, Ar ~ lOkpc (Smith 
et al. 2002), and therefore 

th _ 03 

For the heating time of order the sound crossing time, and 
both much less than the age of the source, the swept-up 
layer of gas should be both isothermal and isobaric. How- 
ever, if the thermal conductivity is significantly suppressed 
over the canonical form given by Cowie & McKee then the 
gas will evolve adiabatically behind the shock. Recently, Et- 
tori & Fabian (2000) have argued, based on the analysis of 
the X-ray data for Abell 2142, that the thermal conductiv- 
ity is very significantly suppressed giving r\ ~ 1/250 to as 
low as ~ 1/2500 at least in some regions of the cluster. If 
the suppression of the thermal conductivity is of this order 
in all cases then the gas behind the shock should clearly 
be treated adiabatically. Given the remaining uncertainly in 
the value for the thermal conductivity, I shall discuss both 
the adiabatic and isothermal cases below. The temperature 
of the swept-up gas will be determined by dynamical pro- 
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Figure 1. The simplified geometry assumed for models of the 
cocoon expansion 

cesses (e.g. shock heating, adiabatic expansion) and radia- 
tive cooling. I will refer to the former processes as dynamical 
heating/ cooling. 

The mass of swept-up gas affects the dynamics of the 
radio source in two ways: firstly the inertia of the gas must 
be considered and secondly the gravitational force on the 
gas due to the cluster potential. As I will show below the 
inertia of the gas does not alter the dynamics in the sense 
that the radio source may expand in a self-similar fashion 
as shown by Kaiser & Alexander (1997; hereafter KA97) - 
their model implicitly allows for the inertia of the gas. A 
simple calculation demonstrates that the gravitational force 
of the swept up gas will also be important towards the end of 
the life time of the source. As the source expansion slows the 
cocoon pressure drops towards pressure equilibrium with the 
confining atmosphere. For gas in a dark matter dominated 
cluster with no pre-heating of the ICM, the pressure in the 
gas is simply determined by the requirement for hydrostatic 
equilibrium in the potential well, hence the pressure in the 
gas must be comparable to the gravitational potential energy 
per unit volume of gas at any point in the cluster. This 
gravitational force must therefore be at least as important as 
the pressure in the atmosphere as the source moves towards 
sub-sonic expansion. 

Throughout this paper I will use the following conven- 
tions: subscripts j, c, h, s will refer to the jet, cocoon, 
hotspot and swept-up gas respectively; the parameters in 
the external gas will be denoted by the subscript x. The ex- 
ternal atmosphere is assumed to follow a King profile of the 
form 

Px = (1 + (r/ao) ) 

which is approximated as a piece- wise power-law (Alexander 
2000), p x = pair/ao)' 13 , with j3 ~ for r < a and (5 = 3f3 K 
for r > oq. Finally, I assume adiabatic indices for the ICM 
of 5/3 and for the radio-emitting plasma of the cocoon of 
4/3. 



3 EVOLUTION OF THE SWEPT-UP GAS 

The discussion of Section 2 shows that it is the expansion 
of the cocoon which is most important in determining the 



evolution of the swept-up gas. A simplified model is then 
sufficient to investigate this evolution and derive scaling re- 
lationships and I therefore approximate the expansion of the 
cocoon as a spherical expansion following the treatment by 
HRB. 

The modelled geometry of the cocoon and bow shock 
is shown in Figure 1. The radius of the bow-shock is R 
and that of the cocoon (or equivalently the contact sur- 
face) XR. The dynamics of the cocoon are determined by 
the energy input from the jet at a constant rate Qo- As 
shown in KA97 and by HRB for this simplified geometry 
(see also Section 4.2) the bow-shock and cocoon expand self- 
similarly with A a constant and the radius of the bow shock 

3 

given by R — Cao (^j^j ~ ^ > where the characteristic time 

to = {a 5 p /Q) 1/3 . 

Two related parameters which characterise the nature 
of the swept-up gas are: (1) the ratio of the cooling time in 
this gas compared to that of the ICM just ahead of the bow 
shock, and (2) the ratio of volume emissivities due to ther- 
mal Bremsstrahlung for the swept-up gas and un-shocked 
ICM. The former depends on the bolometric luminosity and 
therefore, for gas of order a few 10 7 K, scales as T 1//2 /p, 
whereas the appropriate emissivity for a relatively narrow 
band detector is the specific emissivity and therefore scales 



approximately as p T 



-1/2 



3.1 The isothermal case 

I consider first the case when the gas between the bow- 
shock and contact surface is taken to be isothermal; this 
assumption was used by HRB who then calculated the X-ray 
emissivity for a spherical source expanding in a King pro- 
file. Since the pressure is approximately constant through 
the swept-up gas (Section 2) the density, p s , is uniform 
and can be found simply from the total swept-up mass, 
A-kR 3-13 pqclq /(3— /3), and the volume between bow shock and 
contact surface, 47r(l — A 3 )_R 3 /3. The pressure in the swept- 
up gas follows from the jump conditions at the shock which 



I shall assume to be a strong shock, p 3 



where a dot denotes differentiation with respect to time, 
Px = poiR/ao) -13 is the density just ahead of the shock at 
a radius R, and Mo is the Mach number of the bow-shock 
relative to the ICM for a source of radius R. The tempera- 
ture in the swept-up gas, T s , follows from the ideal gas law, 
p — pkT/fi, where \i is the mass per particle: 



Ts 



-. — -TiVlQ ±x 

kp s 4 p s 



The cooling time of both the swept-up gas and the ICM 
just ahead of the shock are proportional to T 1 / 2 /p, hence the 
ratio of the cooling time in the swept-up gas, t cs to that of 
the ICM just ahead of the shock, t cx is 



tcs _ T} /2 p x _ V]**, ( P^V /2 
tcx ps T x /2 2 ° \Ps J 

The density in the swept-up gas is given by p s 

3 ' 
(3-/3X1-A3)' 



Tp x , hence: 



^ = ^((3-/3)(l-A 3 )) 3/2 M 
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Figure 2. The radial dependence of the temperature and volume 
emissivity in the region between the contact surface and the bow 
shock. The two panels show the results for two different atmo- 
spheres (a) /3 = 0.5 and (b) j3 = 1.5. In both (a) and (b) the lines 
have the following meaning: solid line T/T x adiabatic; dotted line 
T/T x isothermal; dashed line e v /e vx adiabatic; dashed-dot line 
isothermal 



(1 - A 3 ) is equal to 15/4(11 - f3) (Section 4.2), hence 
I"*- > 0.05M since f3 < 2. Except for high Mach num- 
bers the cooling time in the swept-up gas is less than that 
in the external atmosphere ahead of the shock, and this will 
certainly be the case as the source slows to near sub-sonic 
expansion. 

Similarly, the ratio of specific volume emissivities is 
given by: 



€vs 
( vx 



T 

± s 



' 1/2 P 2 S 



T~ 1/2 n 2 



V5 V5(3-/3)(l-A3) 



5/2 



Ma 



and — < 125M . These results are illustrated in Figure 
2. 

The enhanced emissivity of the swept-up gas over the 
ICM ahead of the shock leads to an X-ray bright shocked 
shell as discussed by HRB. The mean isothermal tempera- 
ture in the shocked gas differs slightly from the post-shock 
temperature by a factor 5(3 — /3) /(ll — /3) in this approxima- 
tion, however even in a steep atmosphere the factor is close 
to unity, hence the swept-up gas never cools due to dynam- 



ical cooling significantly below the temperature of the gas 
ahead of the shock. In the later stages of expansion as the 
Mach number tends to unity, the cooling time can become 
short compared to that in the cluster gas - if the timescale 
is comparable to the age of the source then significant radia- 
tive cooling of the swept-up gas can occur. Taking a source 
age of 10 7 yrs and a steep atmosphere, then the cooling time 
in the cluster gas would have to be less than of order a few 
10 8 yrs for the swept-up gas to suffer significant radiative 
losses. Such conditions are only possible in the centre of the 
strongest cooling flow clusters. For sources still within the 
core radius of the cluster the cooling time in the ICM would 
have to be nearly an order of magnitude less. It seems likely 
therefore that if the gas is isothermal radiative cooling of 
the swept-up gas is negligible over the lifetime of the source. 

3.2 The adiabatic case 

The situation is very different if the thermal conductivity 
is sufficiently suppressed so that the swept-up gas evolves 
adiabatically. Consider gas which passed through the shock 
when the bow-shock was of radius Ry, where R is the ra- 
dius of the bow shock at the current epoch and < y < 1. 
The immediate post-shock conditions result in a density and 
temperature (p y , T y ) which are simply related by the jump 
conditions to the density and temperature in the ICM; the 
density across the shock rises by a factor of 4, and the tem- 
perature increases by a factor of jgMy, where M v is the 
Mach number of the shock at a radius Ry. As the source 
evolves to a size R this gas must expand and will cool adia- 
batically. The density and temperature of this gas when the 
source reaches a size R (p s (y), T a (y)) are given by the usual 
relationships 

1/7 



Ps(y) 


( Ps_ 


Py 


\Py 


Ts{y) 


( E± 


Ty 


\Py 



1-1/7 



where — is the ratio of the post shock pressure for cocoon 

Py 

radii of R and Ry. In the self-similar model of KA97 (see 
HRB; and Section 4.2) this pressure ratio has a power-law 
dependence on the source size of 

p s ±+£ 

— = v 3 • 

Py 

Assuming an isothermal atmosphere, the Mach number at 
radius Ry is related in the self-similar model to the Mach 
number at radius R by M y = Moj/' 3-2 ''' 3 . Using these re- 
sults it follows that the density and temperature at the 
present epoch for gas swept up when the source was of size 
Ry are given in terms of the external conditions at a radius 
R as follows: 

4(1-/3) 

Ps(y)=4:p x y 5 

, , 5 9 4(3-1) 

T s(y) = J£ T x M%y—s-. 

The distribution of density and temperature between the 
contact surface and bow shock at the present epoch follow 
from conservation of mass. The mass swept up between radii 
Ry to R(y + dy) must now lie between radii Rx to R(x + dx) 
such that: 
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3 p s (y) = 47T1/ 5 ! dyR 3 p x y~ 

Substituting the above results for the density gives a simple 
differential equation relating x and y which has the solution 



15i/ 



(H-/3)/5 



4(11-/3) 

which is valid for A < x < 1 and A 3 = 1 - 15/4(11 - /?). 
These solutions are illustrated in Figure 2. 

These results are in agreement with the numerical cal- 
culations of HA99 and the adiabatic wind models developed 
by Dyson, Falle & Perry (1980). It follows that the cooling 
time of this gas compared to the cooling time of the gas 
immediately ahead of the shock at radius R is given by 

hM = ^ Mo y^. 
t C x 16 

For (5 < 1 the cooling time of the post-shock gas will, for 
even very modest Mach numbers, exceed that of the gas at 
a similar radius in the ICM. For (3 > 1 some fraction of the 
gas will have a cooling time less than the ICM at that same 
radius. For example, taking f3 = 1.5, Mo = 10, all gas swept 
up when y < 0.54 will have a cooling time less than the ICM 
at a radius R, or approximately 40% of the swept-up gas by 
mass. 

The ratio of the specific volume emissivities follows sim- 



64 

= VE V 



ilarly: 

evsjy) = T s (y)- 1/2 Psfa) 2 
e «<" ' I , ' 2 pi 

Unlike the isothermal case these results depend critically on 
the form of the atmosphere into which the source is expand- 
ing and when the gas was swept up. Immediately behind the 
bow shock the X-ray emissivity will be enhanced (unless the 
expansion is highly supersonic). Between the bow-shock and 
the contact surface the emissivity will fall for j3 < 1, whereas 
for j3 > 1 the emissivity rises towards the contact surface, 
the temperature falls and the cooling time decreases. For 
a steep atmosphere we therefore expect to find a layer of 
cooler, dense X-ray luminous gas adjacent to the contact 
surface. 



4 EFFECTS OF THE SWEPT-UP GAS ON 
RADIO SOURCE DYNAMICS 

I have already discussed in Section 2 how the gravitational 
force of the swept-up gas must become comparable to the 
cocoon pressure as the expansion of the radio source slows 
to being only mildly supersonic. In this section I consider in 
more detail the evolution of the radio source itself. There are 
three aspects that must be considered. As first pointed out 
by Scheuer (1974), the contact surface is liable to Rayleigh 
Taylor instability unless the deceleration of the radio source 
is sufficient to stabilise this interface. I shall argue in Section 
4.1 that where observations suggest that the X-ray emitting 
gas is excluded from the cocoon this stability requirement 
may place a useful additional constraint on the dynamics of 
the radio source. Even before the source approaches pressure 
equilibrium the swept-up gas modifies the dynamics of the 
radio source principally by leading to a departure from self- 
similar evolution since the swept-up gas accumulates prefer- 
entially around the cocoon and not the hotspot. In Section 



4.2 I develop a simple dynamical model for this phase of the 
radio source evolution. Finally in Section 4.3 I consider the 
later stages of evolution when the source expansion becomes 
only mildly supersonic. 



4.1 Rayleigh- Taylor stability of the contact 
surface 

The contact surface between cocoon and the swept-up gas 
is liable to become Rayleigh- Taylor (RT) unstable since the 
density of the swept up material exceeds that within the co- 
coon. The surface is stabilised by the deceleration of the 
contact surface (Scheuer 1974). If the total mass of gas 
and dark matter within a radius r is M(r), and assum- 
ing spherical symmetry for the cluster, then the component 
of the gravitational acceleration perpendicular to the con- 
tact surface at a radius r from the centre of the cluster is 
g± — GM(r) cos(ip)/r 2 , where ip is the angle between ra- 
dius vector and the normal to the cocoon. The deceleration 
of the contact surface cannot be measured directly from ob- 
servations, but can estimated in a (weakly) model dependent 
manner by assuming a scaling relation for the advance of the 
contact surface of the form R c — at s , where R c is an appro- 
priate measure of the size of the cocoon. It follows that the 
acceleration of the contact surface can be written: 

R c = 5 3 (5-1)^. 

For the self-similar model of KA97 the exponent 8 — 3/(5 — 
0), and R c and R c , can both be estimated from observational 
data; for the advance of the hotspot R c is simply the length 
of the source, or for the cocoon it is the radius of the cocoon 
assuming approximately cylindrical geometry. 

Observations with CHANDRA show the cocoons of a 
number of radio sources, in particular Cygnus A, to be de- 
void of X-ray emitting gas. For the contact surface to exclude 
the swept-up gas the surface must remain stable to RT in- 
stabilities otherwise fast-growing short wavelength modes 
would disrupt the contact surface allowing the swept-up 
mass to fill the cocoon cavity. The requirement for stabil- 
ity at any point on the contact surface requires the surface 
to be decelerating and: 

ri 2 

\g x \ < \R c \ = \5 3 (8-l)^-\. 

tic 

Smith et al. (2002) give data for the total mass (dark matter 
and gas) within a radius r in the Cygnus A cluster. Using 
these results the perpendicular component of the gravita- 
tional acceleration at the hotspot of Cygnus A is approxi- 
mately 7.4 x 10~ 10 ms~ 2 , and for a point on the cocoon ap- 
proximately three-quarters along the lobe 1.4 x 10~ 10 ms -2 . 
For self-similar models the magnitude of the deceleration 
decreases as the atmosphere becomes steeper. To obtain an 
upper estimate for the deceleration I therefore adopt a lower 
limit to P of ~ 1.5 (Carilli et al. 1994; Smith et al 2002), 

• 2 

giving R c = 0.09J? C /R c - The hotspot advance speed can 
be estimated from spectral ageing arguments assuming an 
equipartition field strength within the lobes. This gives for 
Cygnus A L 3 ~ 0.05c (Alexander Brown & Scott 1984) and 
since Lj ~ lOOkpc 
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\LA = 6.75 x 10" 



L, 



0.05c/ ViOOkpc 



For self-similar expansion the acceleration of the side of the 
cocoon is related to the hotspot acceleration via f c = j^Lj, 

giving for the cocoon f c ~ 0.25Lj. The lobe speed derived 
from spectral ageing is almost certainly an upper limit. Tak- 
ing into account all of the dynamical constraints Alexander 
& Pooley (1996) estimate Lj ~ 0.005c. At best therefore the 
contact surface at the hotspot of Cygnus A is likely to be 
just RT stable. I will return to this point in Section 5. 

In principle, this method provides a way of obtaining 
a lower-limit to the advance speed of a radio source when 
there is evidence that the contact surface is still RT stable; 
the method relies only upon the existence of good X-ray data 
and a simple, but robust, model for the source expansion. 

4.2 Dynamical model with swept-up gas 

In this section I investigate whether the swept-up mass has 
any effect on the dynamics of the source before the source 
ceases to expand supersonically. KA99 show that there is 
significant flow in the swept up gas driven by the pressure 
gradient between the hotspot and the sides of the cocoon. A 
simplified model for the dynamics is therefore suggested in 
which there is no significant effect at the hotspot and only 
the dynamics of the lobe are affected by the gravitational 
force on the swept-up gas in the cluster potential. In order 
to make progress I therefore adopt a simplified approach in 
which I treat the expansion of the cocoon as purely spher- 
ical; this is the approach taken by HRB and Reynolds & 
Begelman (1997). The expansion of the source in the hotspot 
region will be similar to that discussed by Falle (1991) and 
KA97, and the two parts of the overall solution can be re- 
lated by requiring that the jet is in pressure balance with 
the cocoon as in KA97. 

The governing equations are similar to those given by 
HRB, although I now allow for both the inertia and gravita- 
tional force on the swept-up mass and use a more accurate 
treatment of the jump conditions at the shock. The dynam- 
ical analysis assumes the gas behaves adiabatically through- 
out. As before, the radius of the bow shock is defined as R 
and the radius of the cocoon is written as R c — AR. The 
first equation gives the conservation of energy in the cocoon 
assuming that negligible energy is stored in the hotspot: 



7c — 1 dt 



7Pc-T- = Qq. 

7c — 1 at 



Assuming a uniform pressure in the swept-up layer, energy 
conservation in this gas gives 

1 _ . d Ps -y s dV a 15 4tvR 2 ( R \ ~ p - 3 
7s — 1 dt 7 S — 1 dt 16 2 \ ao J 

where the term on the right hand side is just the rate of 
heat input at the bow shock (the factor of 15/16, which 
differs from HRB, allows for the bulk kinetic energy of the 
gas post-shock). The approximation of a uniform pressure 
between bow shock and contact surface is not exact when 
the inertia of the swept-up gas is not negligible. To simplify 
the problem I assume, given the results of Section 3, that 
most of the mass resides in a relatively thin layer either 
at the contact surface when (5 > 1, or just behind the bow 



shock when /3 < 1. Across the rest of the region between bow 
shock and contact surface I assume a uniform pressure of p s . 
The pressure difference between cocoon and the shocked gas 
can be equated to the acceleration of this narrow dense layer 
plus the gravitational acceleration. For j3 > 1 



Ps )=M s ((\R)+g D (\R) 



4ty\ 2 R 2 { Pc 



Where M s = R 3 ^ 13 is the swept-up mass, and <7d(?") 
is the gravitational acceleration towards the centre of the 
cluster at radius r, due mainly to the dark matter in the 
cluster. When j3 < 1 this equation is modified in the sense 
that XR is everywhere replaced by R. This equation differs 
from the equivalent form given by HRB by inclusion of both 
an inertia term and the gravitational acceleration; as out- 
lined below the inclusion of the inertia term still permits 
a self-similar solution to be found when the gravitational 
term is negligible (KA97 implicitly allow for the inertia of 
the swept-up gas in their dynamical solution). 

The gravitational acceleration is found in a self- 
consistent manner by assuming that the distribution of dark 
matter is such as to maintain the gas in hydrostatic equilib- 
rium; for an isothermal gas distribution it follows that the 
gravitational acceleration is given by: 

kT x 1 dp x 

9d = -j-, 

(i p x dr 



and for a King profile 



gD = 3/3k 



kT x r 



aoH a [l + (r/a ) ] 



which can be approximated as follows: 

H if 



9d 



r » ao 
r < a 



where Xu = 3kT x /3pc / pao- 

The final relationships needed to close this set of equa- 
tions are V c = 4ttA 3 _R 3 /3, V s = 4tt(1 - \ a )R 3 /3, and 

p s = |po y^j R 2 > w hich follows from assuming the bow 
shock is a strong shock. In this expressions and throughout 
I will assume 7 S =5/3 and 7 C = 4/3. 

These equations can now be cast in dimensionless form 
by writing t — tor, R — aof(r), where / and r are the 
dimensionless bow shock radius and time respectively. Fur- 
ther, I shall write p c = C{t) Ps and note that given the con- 
stants of the problem (Qo, ao and po) the characteristic time 

to — I -q^- ) ■ In terms of the dimensionless parameters 
the three governing equations become: 

3A 3 / 3 ^(c/ 2 r /3 )+4c/ 2 r /3 ^(A 3 / 3 ) = i 
|(i-A 3 )/ 3 i(f/^) + |/ 2 r^((i-A 3 )/ 3 ) = 



dr 



3(3_^) a2/2 _, /2 = f - ^ (A/) + Yo{XffD 

where Yd = toXo/ao, Pd = ±1 and a dot indicates differen- 
tiation with respect to r. Again the third equation is written 
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for (3 > 1; for < 1 there is no A dependence. These equa- 
tions admit a self-similar solution when the gravitational 
term is negligible Yd « with A and £ equal to constant 
values (Ao, Co) and / = Ct s . Substituting for this form for 
/ gives, after a little algebra, (1 — Ao) = 



4(11-/3) ' 5-t 



and Co 



for /3 > 1 and a similar expression 



S 3(3-/3)A 

omitting the Ao for j3 < 1. The results are in agreement with 
the results of KA97 and the previous section. 

To investigate the case when Yd is small but non-zero, 
I seek a perturbed solution of the form / = /o(r)(l + F), 
A 3 = Ag(l + L) and ( = (a(l + Z), where f (r) = Cr 3 /' 5 "^ 
the self-similar solution. These trial solutions are then sub- 
stituted into the original equations and any terms which are 
quadratic or higher powers of F, L, Z, Yd or any differential 
of one of these terms with respect to r are neglected. The 
remaining equations are then linear combinations of F, L, 
Z, tF, r 2 F, tL and t 2 L, together with the gravitational 
term which has the form t 2 Y£>(Ao/o) /3d_1 . A solution then 
exists such that each term (F, L, Z) has a power-law de- 
pendence on dimensionless time of the form r a , and each 
term has the same time dependence as the perturbation, 



2 Y D (\of ) 



/3d -1 



giving a = 2 + 3(/3 D - l)/(5 - /3). Substi- 



tuting this form for the perturbation into the original equa- 
tions results in a set of 3 coupled linear algebraic equations 
which can be solved to give the magnitudes Li, Z\) of 
each term (e.g. F = Fit") in terms of the perturbation Yd- 
Further details are given in the appendix. For large cluster 
radii, /3d = — 1 and if the atmosphere is steep f3 « 2, a ~ 
and we recover a self-similar solution since the perturbation 
does not grow with time. 

These solutions can be used to investigate how the as- 
pect ratio of the source, which I define here to be given 
by the ratio of the jet to the radius of the cocoon, Lj/R c , 
changes with time. The time dependence of R c in terms of 
the perturbed solution is straightforward R c /ao = A/ = 
^o/o(l + F + L/3) to first order. To determine the time 
dependence of Lj in this approximation we return to the 
method employed by KA97 where the jet as it propagates 
through the cocoon is assumed to be confined by the cocoon 
pressure, pj — p c . KA97 show that this implies that the 
hotspot pressure and cocoon pressure must have the same 
time dependence; this result still holds in the current ap- 
proximation. 

Assuming no accumulation of material at the hotspot, 
the advance of the hotspot is simply determined by balanc- 
ing the ram pressure to the hotspot pressure: 

( L A~ p f 2 

p Lj — p h OC p c 

\a J 

where the time dependence of p c must be given by the 
perturbed solution. For an approximately power-law depen- 
dence of Lj on time, Lj ~ Lj/it^r). The cocoon pressure 
depends on £/~' 3 f 2 , hence writing the time dependence to 
first order 



L (2-g)5 x r (2-p)s-2 



{2-f3)F+^F + Z 



and R c oc r (1 + F + L/3). In Figure 3 I show how Lj and 
Rc evolve with time under the perturbation. The cocoon ra- 
dius changes slowly compared to the length of the jet; the 
cocoon pressure in the perturbed solution exceeds that in 
the self-similar case at the same time since there must now 




10 20 30 40 50 60 70 80 90 100 
t/t Q 

Figure 3. The evolution of the jet length and cocoon radius for 
the perturbed solution discussed in the text compared to the un- 
perturbed solution. The curves are normalized to the unperturbed 
solution at t = to- 
be a greater pressure difference between shocked gas and the 
cocoon to overcome the gravitational force on the swept-up 
gas and the cocoon volume is correspondingly smaller. This 
increased lobe pressure acts to collimate the jet since the jet 
is assumed to be in pressure balance with the cocoon mate- 
rial and the jet therefore propagates more quickly through 
the external atmosphere. The evolution of the aspect ratio 
of the source from its unperturbed (similarity) value follows 
from the above considerations and is given by: 



Rc 



Rc 



L 2aF Z 
~ 3" + 5(2 - 0) + 2-/3 



The perturbation term is proportional to Ydt 01 , and the 
associated coefficient is found by solving the algebraic equa- 
tions as given in the appendix. I now consider these results 
for a typical case when [3k = 0.5, and in the two limits when 
Lj <C ao i-e. (3 — and [3d = 1, and Lj 3> ao i-e. (3 — 3/2 
and (3d — — 1: 



Rc 



■ 0.2Y d t 2 

■ 2.9Y d t 2/7 



Lj <§C a 
Lj » a 



These results can be cast in a more useful form as 2 follows. 
The perturbation term has the form Yd ~ f kTx l a , which 



is equal to 



where c sx is the (constant) sound speed 



in the ICM. Furthermore, since Li 



5-/3 



can be re-written in terms of Lj to give: 

6/5 



hi 

R, 



\Rc 



1 + 
1 + 



0.04 

IT 1 ' 

0.9 

M' 1 



Lj <§C a 
Lj ao 



In this expression M is the Mach number of the hotspot ad- 
vance at the time t. Although the perturbation grows faster 
for Lj < ao, the effect is never significant since the source 
reaches the core radius at a time t ~ to- For Lj 3> ao the per- 
turbation can be significant especially as the source slows to 
being only mildly supersonic. The perturbation grows most 
slowly for steep atmospheres (as expected), however for the 
example considered here even for (3k = 0.5 the gravitational 
perturbation can significantly change the aspect ratio of the 
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source in the sense that older sources will tend to be longer 
and thinner. Of course, the considerations of the previous 
section show that at some point the contact surface must 
become RT unstable, and hence it is unlikely that in the 
majority of cases very long thin sources can be formed by 
this mechanism alone. 

4.3 The final stages of the cocoon 

The latter stages of development of the cocoon are the most 
difficult to treat with simple analytical arguments. The grav- 
itational force due to the swept-up gas becomes comparable 
to the ram pressure as the source slows to being mildly su- 
personic. 

Of particular interest is the ultimate fate of the swept- 
up gas. One possibility is that as the source slows the onset of 
the RT instability (Section 4.1) leads to this gas falling back 
through the contact surface and refilling the low-density re- 
gion previously occupied by the cocoon. Briigen & Kaiser 
(2001) have performed simulations to model the latter stages 
of evolution of a radio source by considering of a light bubble 
of plasma initially in pressure balance with the gas in the 
cluster. They find very different behaviour depending on the 
initial shape of the cloud. Near spherical clouds are quickly 
destroyed by the RT instability leading to mixing of the co- 
coon and ICM. Elongated clouds can form a mushroom-type 
feature which rises through the potential of the cluster in a 
manner first suggested by Gull and Northover (1973). In the 
wake behind these mushrooms colder material is lifted out 
of the cluster. In their simulations Briigen & Kaiser do not 
include the swept-up gas layer I have discussed in this pa- 
per. For the isothermal case of Section 3 when the source 
reaches pressure balance this layer contains mostly gas with 
a shorter cooling time than the ambient gas in the cluster. 
The buoyant phase studied by Briigen & Kaiser will then be 
very efficient at mixing this gas into the ambient material 
at large cluster radii. 

Further discussion of the later stages of the source are 
beyond the scope of the present paper. We are currently 
undertaking simulations to investigate further this phase of 
evolution. 



5 DISCUSSION 

Given the results of the previous two sections it is now pos- 
sible to compare the model predictions with recent X-ray 
data and to consider the evolution of the radio source and 
the cluster itself. 

5.1 Comparison to X-ray observations 

As discussed in the introduction there exist a number of 
studies of the interaction of radio sources with their envi- 
ronment probed by X-ray data. In this section I will focus 
on those published studies of interactions which have suf- 
ficient resolution and sensitivity to probe directly the gas 
between bow-shock and contact surface. 

The CHANDRA results for the Perseus cluster (Fabian 
et al. 2000; Fabian et al. 2002) show clear evidence of a radio 
source displacing the cluster medium. Holes are observed in 
the X-ray emission spatially coincident with the bright radio 



emission from 3C84. Surrounding these holes are bright X- 
ray rims which are cooler than the surrounding gas. Fabian 
et al. (2002) conclude that the source must be expanding 
subsonically since the heating effect of a strong shock is not 
observed. The results of Section 3 provide an alternative ex- 
planation for these data although the explanation proposed 
by Fabian et al. cannot be ruled out. For supersonic expan- 
sion into a declining atmosphere under adiabatic conditions 
I have shown that a layer of cool bright X-ray gas surrounds 
the cocoon, although immediately behind the shock the gas 
temperature is still increased. The observational data will 
effectively measure an emission weighted temperature and 
therefore the cool bright rims surrounding the holes can be 
interpreted as emission from a layer of cool gas adjacent to 
the contact surface. A similar conclusion applies to the data 
for Abell 2052/3C317 (Blanton et al. 2001) in which cool 
X-ray rims are again observed surrounding holes in the X- 
ray emission spatially coincident with the radio structure. 
In this case optical emission is also observed from the X-ray 
rims (Baum et al. 1988) consistent with a cooling time in 
the swept-up gas which is significantly reduced compared 
to that in the undisturbed ICM again consistent with the 
analysis of Section 3. Both of these comparisons should be 
treated with some caution since the radio structure of both 
3C84 and 3C317 is complex and does not resemble a classical 
FR-II, however the analysis presented here is still applicable 
if the radio structures are surrounded by a bow shock. 

The best test case is that of Cygnus A. In this case the 
recent CHANDRA results presented by Smith et al. (2002) 
permit an estimate of the mass of cool gas which may be 
present. For a hotspot advance speed such as that derived 
from spectral ageing of order 0.04-0. 05c (Alexander et al. 
1984), and a ratio of cocoon radius to length of approxi- 
mately 0.2 - 0.25, the expansion velocity of the cocoon is of 
order 0.01c. Using the data from Smith et al. at a radius 
typical of a point mid-way along the cocoon (~ 70kpc) the 
sound speed in the ICM is of order 0.003c, giving a Mach 
number for the sideways expansion of the cocoon of about 
3. The results of Section 3 suggest that all gas swept up for 
y ~ 0.2 will now be cooler than the gas just ahead of the 
bow-shock at a cluster radius of ~ 70kpc. Estimating the 
mass of gas this corresponds to is very uncertain. The linear 
extent of Cygnus A in the cluster is of order lOOkpc and the 
total gas mass within this radius is or order 1O 12 M0; within 
this radius Cygnus occupies about 6% of the volume giving 
a total mass of swept up gas of order 6 x 10 10 Mq. Again, 
assuming a relatively steep atmosphere this suggests that 
approximately 1O 1O M0 of gas swept up by the radio source 
may have cooled to below the temperature of the ICM. As 
discussed by KA99, flow between the bow-shock and contact 
surface will lead to accumulation of gas around the edge of 
the cocoon nearest to the host galaxy and it is here that the 
cool gas should therefore be found. The detailed calculation 
of KA99 when applied to Cygnus A also predicted an annu- 
lus of cold gas in this location. Interestingly, the observations 
of Smith et al. (2002) show conclusive evidence for arcs of 
cooler gas arranged in a belts around the host galaxy which 
are naturally explained as resulting from the accumulation 
of adiabatically cooled swept-up ICM. If the advance speed 
of the hotspot and cocoon are significantly less than this 
value as suggested by considering the ram pressure within 
the hotspot and cocoon (Alexander & Pooley 1996) then 
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the discussion of Section 4.1 shows that the contact surface 
may be RT unstable. If this is the case then the swept up 
gas may well be being entrained within the cocoon. This 
is most likely to occur in the first instance to the sides of 
the cocoon where most cold material has accumulated. In 
this case the belts of cooler X-ray gas will be falling through 
the contact surface and may, as suggested by by Smith et 
al. (2002) provide a mechanism for continued fuelling of the 
AGN. 

5.2 Evolution of a radio source 

The radio source sweeps up a very significant mass of gas 
during its lifetime which is displaced from the centre of the 
cluster. Depending upon the effectiveness of thermal con- 
duction this gas may cool adiabatically to form a thin dense 
layer of relatively short cooling time (relative to the ambient 
ICM at a similar distance from the cluster centre). 

Both the inertia and gravitational forces exerted by the 
swept-up gas will modify the dynamics of the radio source. 
Using the results above it is possible to determine in outline 
the likely evolution of a powerful source. The initial stage 
of evolution will be as the source expands through the core 
of the cluster or the ISM of the host galaxy. In either case 
the mean density of the atmosphere will be approximately 
constant. During this phase the source will follow a self- 
similar form of evolution as discussed by KA97 and Alexan- 
der (2000). Although the effects of the swept-up mass are 
growing, they do not become significant while the source is 
in this phase (Section 4.2). As shown by KA99 there will be 
flow of the swept-up gas away from the region ahead of the 
hotspot to the side of the source. Despite this preferential 
accumulation of material at the surface of the cocoon self- 
similar expansion will still occur since there is a feedback 
between the cocoon pressure and hotspot pressure brought 
about by the jet reaching pressure balance with the cocoon. 

If the source is sufficiently powerful, and the jet is 
not disrupted by a Kelvin-Helmhotz instability (Alexander 
2000), the source will reach a radius equal to the core ra- 
dius of the atmosphere. The hotspot clearly reaches the re- 
gion of declining density first, and as the source passes from 
the cluster core to a declining atmosphere the aspect ratio 
is likely to increase in a phase of non self-similar expan- 
sion. The cocoon pressure will remain high determined by 
the properties of the atmosphere in the cluster core, the 
jet pressure will therefore also remain high and hence the 
hotspot pressure will exceed that predicted for self-similar 
expansion into a declining atmosphere. This phase will cease 
when the bow-shock fully enters the declining atmosphere, 
and the source re-enters a phase of self-similar expansion. 
Beyond this point the effect of the swept-up gas (Section 
4.2) becomes important and as the Mach number of the jet 
decreases the aspect ratio of the source will again increase. 

Throughout this expansion the contact surface may be 
liable to RT instability. This is most severe when the decel- 
eration of the contact surface is small. During the transition 
of the source through the cluster core this is likely to be the 
case. Any instability will allow swept-up gas to be entrained 
into the cocoon material and since the density of this gas 
increases towards the centre of the source, it seems likely 
that any instability will be most apparent in this region. 
This process is also seen in the simulations of Reynolds et 



al. (2001). The source may well survive this entrainment if 
the jet is not disrupted, but this does form a natural mech- 
anism for excluding the cocoon from the region around the 
host galaxy. 

The final stages of the evolution of the source will be 
dominated by the onset of the RT instability as the source 
slows. The gravitational forces on the swept-up gas will be 
largest near the host galaxy (where most of the swept-up 
gas has accumulated). These forces will lead to a pinching 
off of the cocoon near the host galaxy and the onset of a 
buoyant phase (Briigen & Kaiser 2001). It may well be that 
the ingress of material at this time disrupts the jet and effec- 
tively kills the source. The buoyant phase will result in the 
mixing of the remaining swept-up gas as discussed above or, 
when the geometry is optimal, the formation of a true mush- 
room cloud of radio-emitting plasma as discussed by Briigen 
& Kaiser (2001). An alternative scenario is that the RT in- 
stability is effective over the entire contact surface; in this 
case mixing of the radio emitting plasma and swept-up gas 
will be very efficient and the cocoon will be filled with X-ray 
emitting material with the radio-emitting plasma having a 
non-unity filling factor. Pockets of low-density cocoon ma- 
terial will remain buoyant leading to very efficient mixing. 

The analysis of the previous two sections has considered 
a simple model for the cluster atmosphere and in particular 
the relationship between the swept-up mass and the depth 
of the potential well. If the cluster gas were pre-heated the 
gas pressure in the atmosphere will dominate over gravita- 
tional forces of the swept-up layer and the mass of this layer 
may be negligible. In principle it may be possible to use the 
dynamics of radio sources to provide an additional probe of 
the cluster environment, however this requires an indepen- 
dent method of measuring the pressure within the cocoon 
and expansion speed of the source. The latter may be acces- 
sible to VLBI in the future, however at present the equation 
of state of the radio-emitting plasma is very uncertain (e.g. 
Hardcastel & Worrall 2000) and this must be resolved be- 
fore radio sources can be used as detailed probes of their 
environment. 

5.3 The effect on the ICM 

It is now possible to consider the effect of the radio source on 
the ICM. Whether the thermal conductivity is suppressed or 
not, the radio source must do significant work on the ICM; 
the total pdV work must be comparable to the stored energy 
within the lobes (KA97, KA99), and this is deposited as heat 
within the gas between the bow-shock and contact surface. 
However, the effect of the radio source is also to move cluster 
gas out of the centre of the cluster to large radii. 

In the early stages of evolution of the radio source we 
expect high Mach-number expansion into a relatively flat at- 
mosphere. Whether the thermal conductivity is suppressed 
or not the results derived above show that the swept-up gas 
is efficiently heated, a bright X-ray shell develops between 
the cocoon and bow-shock and this gas has a cooling time 
longer than the surrounding material. Furthermore, the ra- 
dio source is lifting a large mass of gas out of the centre of 
the cluster. Clearly during this phase the radio source halts 
any cooling flow in the cluster. 

As the source expands into a declining atmosphere this 
picture changes. If the gas behaves adiabatically and the 
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cluster atmosphere is relatively steep, j3 > 1, then the effect 
of the radio source is to deposit a substantial amount of gas 
at large cluster radii which has a cooling time less than the 
existing ICM at a similar radius. Just ahead of the contact 
surface there will a layer of cooler gas. In the model devel- 
oped above there would apparently always be some gas at 
arbitrarily low temperature, however this is clearly an arte- 
fact of assuming a power-law atmosphere extending to the 
centre of the cluster. In reality, this layer will be thin, while 
the source is expanding highly supersonically, however as 
the Mach number of the expanding cocoon falls the mass of 
cooler gas will become significant. Even for what we think 
of as powerful supersonically expanding sources a significant 
mass of gas may be in a cooler phase. 

As the source slows towards sub-sonic expansion the 
onset of RT instabilities as discussed above will lead to a 
mixing of the swept-up gas with the ambient gas and the 
cocoon material. The effect of this process on the cluster 
may be dramatic. If the swept-up gas has evolved adiabat- 
ically and the expansion has been into a steep atmosphere 
then the effect is to deposit gas at large cluster radii with a 
cooling time less than the ambient material. Indeed some gas 
adjacent to the contact surface may be quite cool and have 
a very short cooling time. It seems likely therefore that the 
mixing process driven by the buoyancy of the radio plasma 
will lead to a multi-phase ICM. As the mixing proceeds this 
multi-phase gas must refill the cavity occupied by the radio 
source. This will either be dominated by fluid-dynamical 
processes or dynamical infall of the densest gas. In either 
case the timescale will be of order the sound crossing time, 
t s , across the radius of the cluster out to which the source 
expanded. The age of the radio source however is of order 
t a /M where M is an average Mach number for the expan- 
sion. Therefore the timescale for the cluster gas to refill the 
cocoon must be of order M times the age of the source when 
it dies; this can clearly be as long as an order of magnitude. 
During this phase their will be no distinct radio source (a 
cluster halo source may possibly exist), but there will be a 
very significant inflow of multi-phase cluster gas. At the end 
of this period the cluster within at least the radius reached 
by the source will be filled with a multi-phase medium with 
a cooling time in the outer parts of this region reduced over 
that prior to the existence of the radio source. 

Similar conclusions concerning the effect of AGN heat- 
ing were reached by Binney & Tabor (1995); their analysis 
also predicts suppression of the cooling flow followed by a 
period of catastrophic cooling after the period of AGN heat- 
ing. The long term effect of a radio source event on the clus- 
ter is however likely to be determined by the buoyant rise 
of material through the cluster atmosphere (Churazov et 
al. 2001; Bruggen & Kaiser 2001; Quillis, Bower & Balogh 
2001). The subsequent evolution of this (no doubt turbu- 
lent) multi-phase medium is difficult to determine without 
the aid of simulations. Whether a subsequent cooling flow 
is enhanced or suppressed will be determined not only by 
the details of this multi-phase ICM, but also by the epoch 
at which the cluster is observed relative to the last radio- 
source event. 



6 CONCLUSIONS 

The effect of a radio source on the cluster gas and vice- 
versa has been considered. I have shown that, if the thermal 
conductivity in the gas is suppressed as has been suggested 
by recent observations with CHANDRA, then a cool dense 
layer of swept-up ICM material will form around the con- 
tact surface of a powerful FR-II radio source expanding into 
a steeply falling atmosphere. This results simply from the 
adiabatic expansion of dense gas as it is carried out of the 
potential well of the cluster ahead of the advancing cocoon 
of the radio source and only occurs when the density of 
the cluster atmosphere falls faster than 1/r for large cluster 
radii. As the source slows from highly supersonic expansion 
to near sub-sonic expansion almost all of the swept-up gas 
has a cooling time less than the ambient ICM; some gas is 
likely to be very cool and there is the possibility of promot- 
ing the formation of a multi-phase medium. 

Eventually the layer of swept-up gas must become 
Rayleigh- Taylor unstable. The swept-up gas accumulates 
preferentially at the edges of the cocoon away from the 
hotspot due to flow in the swept-up layer driven by the 
pressure gradient which must exist due to the difference in 
pressure between hotspot and cocoon. Once the instability 
sets in gas will fall through the contact surface and begin to 
fill the cavity formed by the cocoon. Using existing simula- 
tions as a clue to the behaviour of the source in this phase 
it would seem that depending on the detailed geometry of 
the source at this time there will be both infall of material 
and mixing of the swept-up gas with the ambient ICM. 

The long-term effect of the radio source on the cluster 
when there is (a) suppression of the thermal conductivity 
and (b) efficient mixing of the swept-up layer with the am- 
bient gas, may well be to reduce the cooling time of the 
gas at large cluster radii. This may promote the formation 
of a multi-phase ICM and/or for a period of order a few 
times the lifetime of the source enhance or instigate a cool- 
ing flow. Certainly there must be a substantial flow back 
into the cluster core as the swept-up mass falls through the 
contact surface and fills the cocoon. 

The swept-up ICM must also have a significant effect 
on the dynamics of the radio source especially as the source 
slows towards sub-sonic expansion; if there is no pre-heating 
of the cluster gas then the gravitational force of the swept-up 
gas per unit area of the cocoon must be at least comparable 
to the external pressure. For those sources which escape the 
core radius of the cluster the gravitational force due to the 
swept-up gas can lead to an increase in the aspect ratio 
of the source as the Mach number of expansion decreases. 
For very steep atmospheres falling nearly as fast as 1/r 2 
this perturbation does not grow and the source can expand 
throughout its life in a self-similar fashion. 
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APPENDIX A: THE DYNAMICAL MODEL 

In this appendix I give more details of the solution to the 
dynamical equations governing the evolution of the cocoon. 
The three governing equations in dimensionless form are (see 
Section 4.2): 



d 



d 



sav^ (cfr 



+ 4c./ 2 .r /3 -f (a 3 / 3 ) = - 



(A2) 



3(3-0) a2/2 _, /2 = /3 _, ^ (A/) + Y ^Xf) PD ) (A3) 

where Yd = t^Xo/a® (3d = ±1 and a dot indicates differ- 
entiation with respect to r. The third equation as written 
assumes /3 > 1; for (3 < 1 the equivalent equation is found 
by replacing each occurrence of A/ by /. 

These equations admit a self-similar solution when the 
gravitational term is negligible Yd ~ with A and C equal 
to constant values (Ao, Co) an d / = Cr s . Substituting for 
this form for / in Al gives 

1 27 
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and from A2: 
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4(11-/3) 
and from A3: 

'l-5\ 4 



Co = 1 - 



3(3-/3)A 



and S — 3/(5-/3). This result again holds for /3 > 1, for 
0<1: 



Co = l 



1-5 



3(3 - /3) 



When Yd is small but non-zero, I seek a perturbed so- 
lution of the form / = /o(r)(l + F), A 3 = A u (l + L) and 
C = Co(l + Z), where / (r) = Cr 3/(5 ~' 3) . Substituting these 
trial solutions into Al to A3 and neglecting any terms which 
are quadratic or higher powers of F, L, Z, F, F, L, L, or 
Y D gives: 

= 3[(2 - (3)8 - 2] [L + Z + (5 - (3)F + ^-Fj + 
3r [Z + (2 - (3)F +^F+ ^-F^j + 
125 [Z + L + (5 - (3)F + y F j + 



4r L + 3F 



(Al) 



(A4) 



f*((5-/3)^) = 

^A 3 [(2 - (3)5 + 2] ((5 -(3)F-l + + 
^A. ((2 -flF + + 
f A 3 ((5-/3)F-; + ^F) + 

yA 3 (3F-i) (A5) 
where A3 = 1 — Ao and / = L/ A3 . Equation A3 becomes 
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2rS (*"+§)+ r 2 + f ) + K>(Ao/b) fe_1 (A6) 

for /3 > 1, and where / = Caor 3/(5_/3) . For /3 < 1 the terms 
involving L, L and L do not appear and Ao/o becomes simply 
/o. A solution then exists such that each term (F, L, Z) has a 
power- law dependence on dimensionless time of the form r Q , 
and this time dependence is the same as the perturbation, 
t 2 Fd(Ao/ )' 3d ~ 1 , giving a = 2 + 3(01, - l)/(5 - 0). 

Substituting this form of the solution with F = F\T a 
etc. we obtain a set of algebraic equations for the amplitudes 
F\, Li and_Ki. 

[(18 - 3/3)<5 - 6 + 3a] ( Zi + L x + (5 - + ^Fi 

+ a(Li+3Fi)=0 (A7) 

F! = 2a(3Fi-Zi) + 
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3[a + (7 - P)6 - 2] ^(5 - 0)F 1 -h + ^-F^j 



(A8) 



[2aS + a(a - 1)] (Vi + + yDA^-'c" - 1 (A9) 

Again equation A9 is written for f3 > 1 and the equivalent 
form for /3 < 1 has no terms in Li and Aq 0-1 is replaced 
by unity. Given values for /3 and /3d, the solution of these 
equations is elementary. 



